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In the present work the conditions appearing in the so-called WKB approximation formalism of 
quantum mechanics are analyzed. It is shown that, in general, a careful definition of an approxi- 
mation method requires the introduction of two length parameters, one of them always considered 
in the text books on quantum mechanics, whereas the second one is usually neglected. Afterwards 
we define a particular family of potentials and prove, resorting to the aforementioned length para- 
meters, that we may find an energy which is a lower bound to the ground energy of the system. 
The idea is applied to the case of a harmonic oscillator and also to a particle freely falling in a 
homogeneous gravitational field, and in both cases the consistency of our method is corroborated. 
This approach, together with the so-called Rayleigh-Ritz formalism, allows us to define an energy 
interval in which the ground energy of any potential, belonging to our family, must lie. 



I. INTRODUCTION 



In quantum mechanics the situation is, in a certain sense, quite similar to the situation in classical mechanics, 
namely, the number of physical systems of interest for which the corresponding motion equations can be solved 
exactly is not very large. Therefore approximation methods have been developed and play an important role in 
quantum mechanics. Among these approaches we may mention the WKB, Rayleigh-Ritz, or perturbation methods 
These different ideas have realms of applicability that not always intersect. For instance, WKB is applicable 
to states of the corresponding system characterized by very large values of a certain quantum number 2] , and the 
variational method, also known as Rayleigh-Ritz, is employed to find bounds to the ground state energy of a quantum 
system [2|. 

In the present paper we will focus on the WKB method and analyze the conditions that appear in the implementation 
of it. It will be shown that, in general, they involve two conditions. Usually only one is considered, the remaining 
one is somehow neglected. Indeed, in all text books on quantum mechanics the fact that the method yields a good 
approximation only for the case in which there are several wavelengths between the corresponding classical turning 
points is very carefully explained, and mathematically implemented [2]. 

Nevertheless, there is an additional approximation that is always done, and implicitly introduced, namely, at a point 
on the potential curve a Taylor expansion for the potential is done, and then only a finite number of terms of this 
expansion are considered (usually two terms P, f3|). Then this approximate Schrodinger equation is solved, and with 
it an auxiliary function is defined. This auxiliary function is employed to find the relation between the coefficients of 
the two involved WKB functions, one function on the left-hand side and a second one on the right-hand side of the 
classical turning point. 

It has to be underlined that when this is done, in order to implement the method the solution (stemming from 
the approximate Schrodinger equation) is compared against the WKB wavefunction. This procedure is carried out 
on both sides of the corresponding classical turning point. We must distinguish, in connection with the approximate 
Schrodinger equation and its corresponding solution, that the validity region of the approximate equation (remember 
that it is obtained truncating the Taylor expansion for the potential) does define the region in which the corresponding 
solution can be applied. Of course, the solution to the approximate Schrodinger equation may have a definition 
domain larger than the region in which the approximate potential is valid, but if we use the solution at points outside 
the validity region of the approximate Schrodinger equation, then we use the solution to a potential that does not 
represent very good the physical situation. In other words, a rigorous procedure implies that the solution to the 
approximate Schrodinger equation has to be compared against the WKB wavefunction within the validity region of 
the aforementioned differential equation. 
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This last remark implies that there are two different length parameters involved in the WKB procedure: (i) the 
first one tells us that we cannot be very close to the classical turning point, otherwise the approximate wavefunction 
diverges P, [2] ; (ii) the second length parameter appears when we recognize that the comparison between the solution 
to the approximate Schrodinger equation and the WKB wavefunction has to be carried out within the validity region 
of the approximate Schrodinger equation. Usually, only the first length parameter is considered in the text books 
on quantum mechanics [H, the second one is always neglected, and we may wonder if a careful analysis of it is 
necessary, since we seek an approximate solution, valid at both sides of our classical turning points. 

In the present work we will take into account both conditions and show that a careful analysis of them leads us to 
some interesting facts. For instance, we will show that it is possible to find a lower bound to the ground energy of 
the so-called bound states for some potentials. This is an interesting result, since the current approximation methods 
cannot render a lower bound to the ground state, only the Rayleigh-Ritz method provides a formalism which allows 
us to obtain, via an energy functional, an energy which cannot be smaller than ground energy of the corresponding 
system. With our method the ground energy is always larger than the energy that our method provides. Clearly, 
joining the conclusions of Rayleigh-Ritz with ours we may find an interval in which the ground energy of a one- 
dimensional system has to lie. Though the present approach is restricted to a certain family of potentials we will 
prove that the case of a harmonic oscillator, or a particle freely falling in a homogeneous gravitational field, among 
other potentials, belong to this family. 



II. VALIDITY REGION FOR THE WKB FORMALISM 



The core point in the WKB formalism is related to the fact that the coefficients appearing in the different regions, 
which are separated by the corresponding classical turning points, are matched resorting to the so-called connection 
formulas These formulas are obtained introducing an approximate Schrodinger equation at the classical turning 
point, solving it, finding its asymptotic behavior, and comparing the coefficients against those stemming from the 
WKB formalism. The phrase an approximate Schrodinger equation at the classical turning point plays here a very 
relevant role, and has to be explained in a very clear manner. 

Let us now proceed to do this. Consider now one of the classical turning points, say a;o. A Taylor expansion at this 
point, for the potential is introduced (V^ : / — > 5R, where / C 3?) 

Vix) = Vix„) + V'ix„)ix-x„)+^^^-^{x-xof + .... (1) 

The approximate time-independent Schrodinger equation is defined resorting to the first non-vanishing derivative 
at xq, usually the first-order derivative 



2m dx^ 



+ 



Vixo) + V ixo){x - xo) ^(x) = EiPix). (2) 



The solutions to ^ are Airy functions 1]. The second step is to take the asymptotic form for the corresponding 
solutions to ([2l). This has to be done because the main idea is to use this approximate equation as a patching 
expression, and employ it to connect the expressions stemming from the WKB formalism on both sides of each one of 
the classical turning points. Here the verb connect means to find the relations between the coefficients of the solutions 
obtained with WKB on both sides of the classical turning point. Since this part of the method requires the WKB 
wavefunctions, then the comparison between the Airy functions and the WKB wavefunctions has to be done in a 
region where the latter are well defined, i.e., sufficiently far from the classical turning point [1]. 

Up to this point everything is quite clear in the textbooks on quantum mechanics. Nevertheless, there is an issue 
which is not addressed, and is a relevant one. Namely, the comparison has to be carried out far enough from the 
classical turning point, but resorting to the asymptotic behavior of the Airy functions this last argument is not enough. 
Indeed, we must be sure that in the region where the comparison is done the approximation chosen for the Schrodinger 
equation is a good one, otherwise we would be employing the asymptotic behavior of the solutions of an equation 
which is not a good approximation in the region where the comparison is done. In other words, a careful procedure 
should check that the region in which WKB wavefunctions are valid contains as a subset an interval in which Q is 
valid, i.e., we do not need to take into account the second order derivative, for instance. 

This last statement can be rephrased asserting that the region of the comparison is defined as follows (here from 
the very beginning we assume the usual case, a linear approximation for the Schrodinger equation, additionally recall 
that at the classical turning point E = V{xo)) 
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V{x)^E + V {xo){x-xo). (3) 

If this last approximation is a good one, then the terms of the form {x — a;o)'^, with k > 2, can be neglected and 
this assertion implies 



\V' {x,){x ~ x,)\ > \^^-^ix - x,n (4) 



This last expression defines the validity region of ([2]). Indeed, (|4]) entails that 



In other words, if the matching is done in a region which violates this last condition, then the approximate 
Schrodinger equation should consider higher-order derivatives. Notice that the condition is a function of the in- 
volved classical turning point, xq, in other words, the size of the validity region depends upon xq, it is not, in the 
general case, constant. 

The use of the WKB method requires an additional condition the one can be tracked down to the fact that 
the wavelength of the corresponding particle has to be much smaller than the region in which the potential has a 
noticeable change in its value 

2\E-V{x)\ ^ ^ h 

\dV/dx\ ^2m[E ~ V{x)] ' 

Harking back to ((T|) we impose the condition that 

limsupl y^i)^^^-^ - ^o)| < 1- (7) 

This is tantamount to ask for the absolute convergence of the Taylor series of V{x) Since this last expression 
is valid for all x G /, then limsup] ^ymf^^^J^^ I = 0. This condition is not so stringent as it may seem, for instance, the 
potential of a harmonic oscillator satisfies it, i.e., V'^'^^x) — 0, Vx e 3?. 



III. VALIDITY REGIONS OF THE SEMICLASSICAL APPROXIMATION 



Let us now proceed to define the family of potentials which will be addressed in the present work. 
Firstly, if a; < 0, the potential becomes infinite. Secondly, it is a monotonically increasing function and 



limV{x) 0, as x-^0+ (8) 

Additionally, it possesses bound states, in such a way that the corresponding energy eigenvalues are |ii^07 -E^ii 
where Ei < En, \il < n. 

Let us now suppose that the first and second derivatives of the potential do not vanish Va; G 3?+. This condition 
is introduced in order to have in our family the case of a truncated harmonic oscillator. The case of a particle freely 
falling in a homogeneous gravitational field will also be addressed here, though for this potential the second derivative 
vanishes we may analyze it. The possibility of resorting to the WKB method, using a linear approximation for the 
potential energy at the classical turning point xq , is feasible at those points x such that they fulfill (O and ^ . These 
two conditions imply 
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2.V"{xo) Xo-X 

OL > 



■ 2 a 
Here we have defined the length parameter 



-a- 



2/3 



Axq — X . 



(9) 



a = 2|y'(a;o)/y"(xo)|. 



(10) 



This assertion may be rephrased stating that WKB can not be apphed at a point x in the classical region (the one 
is defined by the interval J = (0, xq]) if the following inequality holds 



2.V"{xo) Xq-X 

a ^ < 

2 a 



V"{xo) 
8m" 2 



2/3 



Axo_ 
3 c 



(11) 



If the use of a linear approximation to the potential is inconsistent we could try to employ the WKB formalism 
considering higher-order derivatives in the approximate Schrodinger equation, and we would end up with an inequality 
(though a more complicated one) and a new length parameter, the one would be given by the expression a = 
(? + l)|y«(xo)/F('+i)(xo)|. 

Let us now analyze the implications of the breakdown of the linear approximation, and study its consequences, if 
any, upon the energy eigenvalues. 
Defining 



{Ama^V'ixo) 



(12) 



We may cast pT|) as follows 



xo — X 
a 



|2| <7|l-4z/3|. 



We consider now the classical region J, and restrict x to it, then z E [0, Xq/q] 
This last condition allows us to define the following two functions 



(13) 



(14) 



3?, 



(15) 



g : [0, So/a] 3?, 



(16) 



(17) 



Kz) = 7|l-4z/3|. 



(18) 



It is readily seen that G [0,xo/a] we have that f{z) > 0. On the other hand, our function g{z) may change sign 
in this interval. Indeed, 



7|l-4z/3| 



7[l-4z/3], when z< 3/4 
7[-l + 4z/3], when z> 3/4 



(19) 



We have two possibilities 
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A. First Case {xo/a < 3/4) 



Then the condition that iniphes the breakdown of WKB reads 

2<7[l-4z/3]. (20) 
The border that divides the regions, in which WKB is vahd from the one in which it is not vahd, is given by 

z = 7[l-42/3], (21) 

1 + 47/3' ^^^^ 

In other words, in the region z E [0, i^J^/^ ] WKB cannot be used, and in z G [ -^^J^^^ jXq/o] WKB is vahd. 
Notice that if the interval z e [ 1+4-^/3 ' ^o/q] becomes one point, then in the region x G (0,a;o] WKB cannot be 
used. This happens when 

-xo/a, (23) 



1 + 47/3 

This condition can be cast as follows 



8m 2 

The roughest approximation allow us to rewrite (j24p as 



f^^v"ix,)\l-'^ . (24) 



3a 



8m 

Additionally, we have the condition xo/a < 3/4, which implies 



^=xlV'{xo). (25) 



2V"{xo) 

Consider now the following potential (/3 > 0) 



^0 < T^TTTTTZr-^- (26) 



V(x) = <! ^^'^^ ^ ^ (27) 

^ ' I 00, when a; < ' 



([^ implies that for this kind of potentials 



n < 5/2. (28) 



Notice that the fulfillment of ([26|) for the aforementioned family of potentials does not involve the value of (3, only of 
n. Under this condition we may find some important potentials of quantum physics, namely, the harmonic oscillator 
{n = 2) and a particle freely falling in a homogeneous gravitational field {n = 1). 

The value of xq, using ((25|) and ([27)l . is given by 



\8pmn/ 
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From this last expression, additionally, we find the following energy 



4'^-/3^o=5 2- (30) 

The meaning of this energy is the following one. For this kind of potential this is the minimum energy that a 
particle can have in order to, with a linear approximation to the potential, be able to employ WKB. Indeed, notice 
that xq is the smallest value of the coordinate at which the aforementioned formalism can be employed, and since we 
have assumed, from the very beginning, a monotonically increasing potential, then we deduce that if a particle has 
an energy smaller than ([3D]), then WKB cannot be employed. 

At this point we may wonder if this energy has some physical meaning. We will show that it is a lower bound 
to the eigenenergies of the corresponding bound states, i.e., all the eigenenergies stemming from the solution to the 
corresponding Schrodinger equation, and related to bound states, are larger or equal to e'^^"' . This assertion will be 
proved in section IV. 



B. Second Case (xo/a > 3/4) 



In this case the breakdown of WKB involves expression (fT9|) . and clearly the first point of interest is related to 
the fulfillment of expression (|2ip . In other words, once again, in the region z E [0, i^J-y/a l WKB cannot be used. 
But, since we have that xo/a > 3/4, now the breakdown of WKB can include a new interval, the one is absent when 
Xo/o: < 3/4. Indeed, consider now the possibility in which 



z = 7[-l + 4z/3]. (31) 

This happens when 

(32) 



-1 + 47/3' 



Clearly, the case 47/3 = 1 has to be discarded in this last equation. Nevertheless, we must also analyze the case 
in which 47/3 = 1. Expression ([20]) tells us that in the region z G [i+^^j^^o/a] WKB can be used, this is when 
47/3 = 1 This is easy to understand since the slope of the two straight lines, f{z) = z and g{z) — 7[— 1 + 4z/3], is 
the same, 47/3 1. 

Harking back to the case 47/3 7^ 1, we conclude that if < xg/a, then we have three different regions, 

namely: 

(i) If z e [0, Yq:^^] WKB cannot be used; and it is related to the fact that we are too close to the classical turning 
point and in consequence ^ is not valid. 

(ii) If ^ e [ 1^27/3 ' 47/3-1 ] WKB can be used; 

(iii) Finally, when z e , xq/o] WKB loses, once again, its validity, but this time the breakdown of the method 
emerges because the approximation introduced for the Schrodinger equation at the classical turning point is violated, 
see expression (ITU)) . 

Of course, this last possibility occurs only when 

This last condition can be cast as follows 



Sm - \ ' J 2 



(34) 
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Consider now the case in which the potential is given by (|27p. One of our previous results, see (l28l) . imposes as 
condition n > 5/2. Then we obtain as condition for the existence of the region defined by z € [ 4.y/3_i ; ^o/o^] the 
following inequality 

/ 2 n3 

7T > 



{^^)nVi.o). (35) 

This last expression allows us to define the following energy 

C-^Y^- (36) 

V 2 / Smnx?, ^ ' 



(1) _ (i\n^ 



C. Comparison between the two cases 



If xo/a > 3/4, then we cannot obtain the value of n without knowing also the value of /3 (the one appears in V{xo)), 
as expression ([55]) clearly shows. This does not happen in the first situation {xo/a < 3/4), see expression in 
which the condition upon n is independent from the value of /3, though xq does involve /3, see (j29p . 

An additional difference between these two cases lies in the fact that if {xo/ct < 3/4) then we may find a family 
of potentials such that the analysis contains only one region, see expression ([24| . in which WKB is not valid (in the 
interval (0, xq) the conditions imposed by WKB cannot be satisfied), whereas in the other situation the analysis shows 
the possibility of having two or three regions, and in one of them WKB can always be employed. 

As a byproduct we have obtained for the first case an energy, see expression ([5(11) such that for energies smaller 
than it the WKB method cannot be used. Similarly, for the second situation, expression p6|) defines an energy such 
that if the energy of our particle is smaller than that provided by ([36]) . then WKB cannot be used. In other words, 
in both cases we find a lower bound for the validity of the method, the one includes a linear approximation for the 
Schroodinger equation at the classical turning point. 



IV. LOWER BOUND FOR THE GROUND STATE ENERGY 



Let us now consider the classical energy E^j^ = V{xo) obtained previously, see expressions ((30| and (|36p. The 
meaning of this energy is the following one. For the kind of potential under consideration this is the minimum energy 
that a particle can have in order to, with a linear approximation to the potential, be able to employ WKB. Indeed, 
notice that xq is the smallest value of the coordinate at which the aforementioned formalism can be employed, and 
since we have assumed, from the very beginning, a monotonically increasing potential, then we deduce that if a particle 
has an energy smaller than (f30|) or ([36]) . then this energy can not be obtained via WKB. 

At this point we may wonder if this energy has some physical meaning. Let us now answer this question. We will 
conjecture that if 

^S' <- 

then it is a lower bound to the eigenenergies of the corresponding bound states, i.e., all the eigenenergies stemming 
from the solution to the corresponding Schrodinger equation, and related to bound states, are larger or equal to E^j\ 
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Let us now consider the lower energy eigenstates of a quantum system subject to a potential satisfying our require- 
ments. Under this condition WKB may be used, though it provides energies that are not a good approximation to 
the correct ones. We now proceed to prove that when E^,]"^ satisfies p7|) . then it does define a lower bound to the 
ground energy of the system. Consider the energy Eq of the ground state and assume that it intersects the potential 
at Xg, i.e., Eq = V{xg). Here we have two possibilities: 

(i) Xg lies on the right-hand side of xq. In this case Eq can be obtained via WKB, and since iJ^j^'' is the smallest 
energy that allows the use of this approach, then £'^|-' < Eq, and we are done; 

(ii) Consider the remaining possibility, Xg < xq. Since we have assumed a monotonically increasing potential, then 
this last condition implies that V{xg) < V{xq). 

The idea here is to proceed by contradiction, namely, assume that Xg < xq, a fact that implies V{xg) < V{xa), and 
resorting to the fundamentals of quantum mechanics we will show that this assumption leads to an inconsistency for 
those potentials fulfilling ((37)) . 

According to quantum mechanics the distance that the wave function (of a bound state) tunnels inside the classically 
forbidden region becomes larger as the energy grows. Indeed, in the aforementioned region the wavefunction behaves 

as i>{x) « exp{-i /; ^2m[Viz) - E]dz'} 0. 

Let us now recall the uncertainty relation 

AxAp>^. (38) 

For the ground state the tunnelling distance is the smallest one, and if Xg denotes the intersection of the ground 
energy with potential, then (0, Xg] is the classical region, and we may state that 

Ax - Xg. (39) 

This last expression says that the probability of finding the particle in the classically forbidden region is almost 
zero, and in consequence, the uncertainty region associated to the position of the particle is related to the classical 
region, in which the wavefunction has an oscillatory behavior. Bearing this in mind the uncertainty relation can be 
cast in the following form 



Xg^2m[Eo- <V>]-{<p >)2) > -. (40) 

Where the averages are calculated with the wavefunction of the ground state. Since we have assumed that V{x) > 
0, Vx e then (gO]) imphes 



Xg^2mEo > -. (41) 
Then the ground state energy satisfies the condition 



a 

And since our initial hypothesis has been the condition Xg < xq, then 

< J}— (43) 
Sma^Q ~ SmXg 

Nevertheless, we have assumed from square one the validity of ([57]) . and hence we find a contradiction. Indeed, our 
three initial conditions (expression (|37p and Xg < xq and the uncertainty relation) imply £'^|'' < Eq. Nevertheless, 
since the potential has been assumed monotonically increasing, then Xg < xq ^ ^i^g) < ^(^^o), and in consequence, 
we find that the simultaneous use of these three conditions is, logically, inconsistent, and therefore the case Xg < xq 

has to be discarded, if expression ([37)) and the uncertainty relation are to be kept. Clearly, we conclude that £'^|^ 
provides a lower bound for the ground energy. We now proceed to analyze if any of our two possible cases satisfies 

dSZl). 
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A. First case 



Consider now a potential given by dail), with 1 < 71 < 5/2. Then, E^.]^ = sS;;^' ^^'^ condition ^ is fulfiUed. 

In other words, if 1 < n < 5/2, then sj^j"^ does provide a lower bound to the ground energy of the corresponding 
system. 



B. Second case 



The remaining possibility, n > 5/2, cannot be analyzed in the context of our approach. Indeed, notice that now we 
require the fulfillment of the following condition 



^ - V 2 J Smnxl ^ ' 

This last condition can be cast as 



\2) - {n-l) 



< 7^^- (45) 



Clearly, not only has to be satisfied, additionally n > 5/2. A fieeting glimpse shows that this is impossible, 
in other words, we cannot use the present approach to find a lower bound for potential of the form V{x) = /3a;", if 
n > 5/2. 



V. SOME PARTICULAR POTENTIALS 



Notice that under the condition n < 5/2 we may find the potential associated to a harmonic oscillator, or a particle 
freely falling in a homogeneous gravitational field. We know proceed to calculate the bounds for these two potentials 
that our method provides. We use these two cases to corroborated our predictions, since we need potentials whose 
exact eigenenergies are already known. 



A. Harmonic Oscillator 



A particular case of the situation just mentioned is the potential associated to a harmonic oscillator. Notice that 
for this case the condition n < 5/2 is fulfilled, and in consequence, if our argument is valid, then we expect to obtain 
a lower bound to the ground energy. 



V(x) = I > ° (46) 

^ ' [ oo, when a; < ^ ' 

In this case (see expression ([^ l 

^0 = ily^'Xo, (47) 
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Ao = \/— , (48) 
V muj 



The eigenenergies of a harmonic oscniator read i?„ = huj[{2n + 1) + 1/2], in this case n is an integer. Clearly, since 
we have the truncated harmonic oscillator potential, only the odd wave functions are allowed and, in consequence, 
the ground state is given by £^o = Here we have that, indeed, E^^' < Eg. 



B. Free Fall 



Another interesting situation in quantum mechanics, the one has exact solution, is the case of a particle falling in 
a homogeneous gravitational field ^5^] . Here the potential satisfies the condition 



Vix) 



mgx, when x > 
oo, when x < 



(50) 



The eigenenergies are given by (n g N) 



En 



fr r37rr li-i^/J 



2mP I 4 



2/3 



(51) 



1 2m'^g 



(52) 



For this type of potential the second derivative vanishes in the region (0, oo), and hence, the condition that deter- 
mines the minimum distance xq is given by 



If fi^ \l/3 
^°-2l^ 



(53) 



Since the potential is linear, then, clearly, the linear approximation to Schrodinger equation is always valid, and 
the breakdown of the method is only related to the violation of ^ . 
The ensuing energy reads 



e[1^ = 



/32m/2 



Notice that the energy eigenvalues can be cast as follows 



En = Ei, 



1) rSTT 
' I 2 



2n 

2 



|2/3 



(54) 



(55) 



The ground state corresponds to rt = 0, and then ([55]) becomes 



(56) 



Clearly, once again, the condition, i?o > E^^^^ is fulfilled, |^^| 



2/3 



> 1. 
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VI. CONCLUSIONS 



In the present work it has been shown that the so-called WKB approximation method involves two length pa- 
rameters which have to be introduced in the calculation of the relevant variables. Usually, only one of these length 
parameters is considered in the text books on quantum mechanics P, 0] , the second one is always neglected. These 
parameters are related to different conditions that have to be fulfilled in order to have a mathematically consistent 
approximation method: (i) the first one tells us that we cannot be very close to the classical turning point, otherwise 
the approximate wavefunction diverges P, 01j see expression ([S]); (ii) whereas the second length parameter appears 
when we recognize that the comparison between the solution to the approximate Schrodinger equation and the WKB 
wavefunction has to be carried out within the validity region of the aforementioned motion equation, see expression 

Both conditions have been taken into account and it has been shown that a careful analysis of them leads us to 
some interesting facts. For instance, it is possible to find a lower bound to the ground energy for some potentials, 
i.e., those whose behavior goes like V{x) = /3a;", /? > and 1 < n < 5/2. Joining the present conclusions with the 
so-called Rayleigh-Ritz method allows us to find an interval in which the ground energy of a one-dimensional system 
has to lie. Additionally, some examples of the approach have been explicitly calculated, for instance, the case of a 
harmonic oscillator, or a particle freely falling in a homogeneous gravitational field, and it has been shown that the 
energy deduced with our arguments is, indeed, smaller than the ground energy of the corresponding system. 
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